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This paper applies the theory of cylindrical magnetic domains 1 to 
cylindrical domain devices. The stability conditions are examined as 
bounds to the region of possible device operation and it is found that 
the plate thickness, h = 41, and the domain diameter, d = 81, where I 
is the ratio of the wall energy per unit area to fa times the saturation 
magnetization squared, are preferred values. When the effects of wall 
coercivity and mobility are examined, it is found that the preferred 
plate thickness and domain diameter are even more strongly preferred, 
that the wall motion coercivity should be less than one percent of ^ 
times the saturation magnetization, and that a domain coercivity and 
mobility may be defined. Consideration of the Neel temperature and 
the desired absolute domain size in addition to the static stability 
conditions shows that domain materials having some antiferromagnetic 
character and induced uniaxial anisotropy are preferred. Where appro- 
priate, domain methods for measuring material parameters are 
described. 

I. INTRODUCTION 

The application of cylindrical magnetic domains or "bubbles" to 
memory and logic devices has recently received considerable atten- 
tion. 2-0 Such domain devices may operate in a continuum of modes 
ranging from the wall motion coercivity dominated mode to the "hard 
bubble" mode. In the coercivity dominated mode, applied fields deter- 
mine the domain configuration which is then maintained by coercivity. 
In the hard bubble mode the coercivity must be sufficiently low that 
the domains have a well-defined size and shape permitting the move- 



* Portions of this article were presented at the "Fifteenth Annual Conference on 
Magnetism and Magnetic Materials" Philadelphia, Pennsylvania, on Nov. 20, 
1969. (See Ref. 2.) 
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ment of individual domains as distinct entities. A previous paper 
developed the theory of static stability of cylindrical domains in 
materials having zero coercivity. 1 The present work applies this theory 
to the determination of the preferred conditions for construction of 
devices operating in the hard bubble mode. It is found that specifying 
an operating domain diameter determines preferred values for the 
plate thickness, magnetization, anisotropy constant and the maximum 
allowable value of the wall motion coercivity. 

Figure 1 shows the model for the domain structure from which the 
static stability theory was developed. The coordinate system and 
symbols used here are the same as in Ref. 1 except for the addition 
of a few symbols such as /x 10 , the wall mobility; H c , the wall motion 
coercivity; v rf , the domain velocity; and Ao-, c the variation in wall 
energy. The model represents a single isolated domain in a plate of 
magnetic material of uniform thickness, h, and an infinite extent in 
the plane, r f = oo. Everywhere within the material the magnetization 
has a uniform (saturation) magnitude, M 8 , directed along the upward 
plate normal (the z direction) within the domain and along the down- 
ward plate normal elsewhere within the material. The domain wall is 
assumed to have negligible width and the domain wall energy density, 
v w , is initially taken to be independent of both wall orientation and 
curvature. The domain wall is cylindrical in the sense that it every- 
where contains a line parallel to the plate normal. Under these assump- 
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Fig. 1 — Domain configuration and coordinate system. 
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tions only the component of the spatially uniform applied field which 
lies along the plate normal will interact with the domain so that only 
this component is considered. This component is denoted by H and is 
taken positive when directed upward, the direction tending to collapse 
the domain. 

The assumptions implicit in this domain model are not all inde- 
pendent. The interrelation of the assumptions and the dependence of 
the assumptions on the domain geometry and the material parameters 
are discussed on pp. 3312-3318 of Ref. 1. The validity of the assump- 
tions will not be discussed further here except to note that in order 
for domains of the type to be considered here to exist, domain nuclea- 
tion and wall-width considerations require roughly that 

K u > 2tM] , (1) 

where K„ is the uniaxial anisotropy constant, and that the validity of 
the approximations generally improves as K„ is increased above this 
minimum value. 

The domain radius function, r b {6) , which is expanded in the series 

r b (6) = r„ + Ar„ + £ Ar n cos [n(0 - 6 n - A0„)], (2) 

n = l 

describes the domain shape in the plane. The Ar„ and A0„ describe a 
variation in domain size and shape from a circular domain of radius 
r b (d) = r and the 6„ describe the direction of the variation. (The A0„ 
have significance only for second variations.) Since only near circular 
domains are of interest here, the condition 

|r„|» |Ar | + f> |Ar„| (3) 

n=l 

is imposed to assure that the radius is single valued and smooth. 

The first and second variations of the total domain energy with 
respect to the Ar„ and A0„ determine the domain equilibrium and 
stability conditions. The total domain energy, 

E T = E w + E„ + E M , (4) 

is the sum of three terms. The wall energy, E w , is the product of the 
wall energy density and the wall area. The applied field interaction 
energy, E„, is proportional to the product of the domain volume and 
the external field interaction energy. The last term, Em, is the internal 
magnetostatic energy of the domain. The energy variation has the 
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where in the energy derivatives the independent variables, the Ar„ 
and A9„ , have been written as r n and 6 n for compactness, the zero sub- 
script indicates that the derivatives are to be evaluated for a strictly 
circular domain, r b (6) = r , and 3 indicates terms of order three or 
higher. The first derivatives are the negative of the generalized forces, 
— (dE T /dr„) and — {dE T /dd„) Q being respectively the radial and angular 
generalized forces of rotational periodicity n. The second derivatives 
form the stiffness matrix of the system with (d 2 E T /dr n dr m ) being the 
(n, m) element of the radial submatrix and (d 2 E T /dr n d6 m ) and 
(d 2 E T /dd n dd m ) being the corresponding elements of the mixed and 
angular submatrices respectively. In Ref. 1 the derivatives of the 
total energy with respect to the A? - „ and A0 n were obtained by differ- 
entiating the integrals which form the terms of equation (4) and evaluat- 
ing the resulting integrals for the case of a strictly circular domain. 
In the present work it is convenient to write the energy variation 
expression in a normalized form in which: energy is measured in units of 
4(27ril/^)(7r/i 3 ), the equilibrium condition has been used to eliminate 
the applied field from the second variations (the stiffness matrix is of 
interest here only for a domain in equilibrium), the first- and second- 
order variations which are identically zero and deleted, and the non- 
zero generalized forces and stiffness matrix elements are written as 
functions of dimensionless variables, the applied field measured in 
units of the magnetization, H/AirM, , and ratios of the plate thickness, h, 
the domain diameter, d = 2r , and the characteristic length of the 
material, 



I = tr w /4t*M\ 



(6) 



This normal form of the energy variation is [Ref. 1, equation (68) 
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where F and the <S„ are called the force and stability functions respec- 
tively and are written in terms of the complete elliptic integral of 
the first kind, 

«r/2 

K(m) = / (1 - m sin 2 0) _i dd, (8) 

■'0 

the complete elliptic integral of the second kind, 

fT/2 

E(m) = (1 - m sin 2 0)* dd (9) 

Jo 

and the L„ functions which are auxiliary functions introduced for 
convenience, 

L (x) = (10a) 

Ux) = 4[(x + l)»fl((l + xT 1 ) - x(x + 1)"*X((1 + xT 1 )] (10b) 
1 



L n +i(x) = 



[4n(2.r + l)L»(x) - (2n - !)£,._, (a:) 



2n + 1 

- 16n.r(x + 1)"* X K((l + a)" 1 )], n fc 1. 

The force and stability functions are 

F(x) = -x\{\ + 0*^((1 + X*f l ) ~ 1], 



S (s) - F(s) -*^F(x), 



= ~x 2 [Li(O-^i(0)], 

27T 



(10c) 

(11) 
(12a) 
(12b) 



and 



1 1 



n(x) = —2-^-7 ^r .T 2 [L n (x' 2 ) - L x {x-*) - LM + L,(0)] f 
n — 1 /7T 



n > 2. 



(13) 



Figure 3 of Ref. 1 plots F(d/h) and S„(d/A) for d/h ^ 6 and n ^ 10. 
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Appendices A and B of Ref. 1 give methods of computation and series 
expansions of these functions while Section IV discusses the physical 
interpretation of the terms in the energy variation expansion [equa- 
tion (5) here]. 

Section II of this paper discusses the domain size and stability 
implications of the energy variation expansion (7). This discussion 
yields the conditions for the existence of cylindrical domains in the 
total absence of dissipative processes and these existence conditions in 
turn place several restrictions upon device design. The section con- 
cludes with a discussion of domain size and domain ellipticity in the 
presence of anisotropic wall energy. Section III considers the effects 
of dissipative processes (wall motion coercivity and mobility) on 
domain existence and movement, the relation of domain mobility to 
wall mobility, and the limiting conditions under which the hard bubble 
mode may be achieved. Section IV combines the results obtained here 
and in Ref. 1 to obtain several relations between material parameters 
and device performance. 

II. CYLINDRICAL DOMAIN SIZE, SHAPE AND STABILITY AT EQUILIBRIUM 

This section examines the domain equilibrium and stability condi- 
tions and some of their implications. 

2.1 The Equilibrium Condition 

The domain is in equilibrium when all of the first-order variations 
of the total energy with respect to the Ar n and A0„ are zero. In the 
variation expansion (7), all the first-order energy variations except 
the variation with respect to Ar„ (representing a variation in domain 
size with no variation in shape) are identically zero. Since the domain 
is initially assumed to be a circular cylinder and there are no forces 
tending to deform it, the domain is in equilibrium when it is a circular 
cylinder having a diameter which is a solution to the normalized force 
equation, 

1 , d _H_ Jd\ . , 1A . 

The (normalized) generalized forces appearing in this equation have 
a one to one correspondence with the terms of the energy sum (4): 
the first term being produced by the wall energy, the second by the 
applied field, and the third by the internal magnetostatic energy. The 
normalized wall force, —l/h = —(T w /4irM 2 ,h always tends to collapse 
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the domain (a w and h are positive) and is independent of domain di- 
ameter. Each of the normalized generalized forces may be converted 
to an equivalent field by multiplying by MicMJd. The equivalent 
wall field is a w /2r M. so that the wall field is proportional to the product 
of the wall energy and the wall curvature. The normalized applied field 
force, -(d/h){H/±irM,), tends to collapse the domain for positive 
applied fields and its magnitude is proportional to the domain diam- 
eter. In this case the equivalent field is the applied field, H. The normal- 
ized internal magnetostatic force is 4- F(d/h). The force function, F, is 
positive with positive first derivative and negative second derivative 
at all points. It approaches d/h for small domain diameters and 
it' 1 In | ie^d/h | for large domain diameters (In Ref. 1 see equation (138) 
and Fig. 3). Since F is positive, the internal magnetostatic force always 
tends to expand the domain. The internal magnetostatic force and the 
wall force are thus oppositely directed. Therefore equilibrium is at- 
tained for any given diameter by adjusting the applied field to a value 
which compensates for the difference in magnitude of these two forces, 
the sign of the field depending on which force is dominant. The equiva- 
lent field for the case of the internal magnetostatic force, 4tM, 
• (h/d)F(d/h), is the z averaged z component of the internal demagnetiz- 
ing field or the internal magnetostatic scalar potential difference be- 
tween the top and the bottom of the plate divided by the plate thick- 
ness. 1 This field approaches 4irM, for small domain diameters and 
4M t (h/d) In | 4e*d//t | for large domain diameters. 

Solutions to the equilibrium problem may be discussed either in 
terms of the equivalent fields as was done by A. H. Bobeck 3 or in 
terms of the generalized force equation (14) with the preferred method 
depending on the specific application. In the present case, the generalized 
force equation will be used since the general properties of the solutions 
of the equilibrium problem may be easily obtained by straight line 
constructions on a plot of F. 

Figure 2 shows examples of such constructions (along with stability 
constructions which will be explained later). The equilibrium con- 
struction consists of first locating the point on the vertical axis whose 
ordinate is l/h, then drawing a straight line through this point whose 
numerical slope is H/4:irM a . The line so constructed thus represents the 
first two terms in equation (14) so that its intersections with F (if any) 
are the equilibrium points. The dashed line asymptotic to the force 
function at the origin has numerical slope one and therefore provides 
a reference for estimating the magnitude of applied fields. In each of the 
constructions of Fig. 2, l/h is 0.3 and this point on the vertical axis is 
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Fig. 2 — Construction of solutions to the force equation for l/h = 0.3, d = do , 
d = cU and d = 2h. 



denoted by W. The symbol V? was used since l/h is the wall energy 
per unit area in units of twice the magnetostatic self-energy per unit 
surface area of an infinite plate of uniform thickness, h, when it is 
uniformly magnetized in a direction perpendicular to the surface of the 
plate. 

For the first example, consider the line WCt' for which l/h = 0.3, 
H/4irM, = 0. This line has one intersection with F at Q(d/h = 0.378). 
Since the slope of F is everywhere positive, there will in general be 
only one intersection with the construction line for any zero or negative 
value of H/^tM, independent of the magnitude of l/h. That this 
single solution is unstable may be appreciated by considering small 
variations in the solution diameter from its value at G: Increasing the 
diameter of the domain from the diameter at 6 increases the magnitude 
of the internal magnetostatic force term while the other terms remain 
constant. The total force tending to expand the domain thus increases 
from zero as the domain expands from the solution diameter. The domain 
is thus unstable with respect to expansion. An entirely similar argu- 
ment shows that the domain is unstable with respect to collapse. It is 
easily seen in general that for a fixed solution diameter the solution 
will become more unstable (the destabilizing force increasing faster with 
increasing diameter) as the field is made more negative. Stable solu- 
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tions to the equilibrium problem exist only for fields aiding the wall 
field in tending to collapse the domain although one unstable solution 
exists for any magnitude of negative applied field. 

Because of the logarithmic behavior of the force function at large 
domain diameters, a radially stable solution to the force equation will 
appear (in addition to the radially unstable solution) when a vanish- 
ingly small positive bias field is applied. (Radial stability is determined 
in the same way as in the preceding paragraph. The existence of radial 
stability does not assure total domain stability as will be seen.) Since 
the F function has a negative second derivative everywhere, there can 
be no solutions in addition to a single radially unstable and a single 
radially stable solution. This situation is illustrated by the line WS 
in Fig. 2. The line represents l/h = 0.3, H/4vM. = 0.254 for which 
the unstable solution, denoted by C U, occurs at d/h = 0.767 and the 
stable solution, denoted by S, occurs at d/h = 2.000. Bobeck 7 has ob- 
served experimentally the unstable solution under static conditions 
which permit the existence of the stable solution. This was done by 
applying a short duration bias field pulse which reduced the domain 
diameter from the stable solution diameter at S to the unstable solu- 
tion diameter at 'U and then returning the bias field to its original value 
just as the unstable solution was attained. The domain having the 
unstable solution diameter was then stabilized sufficiently by coercivity 
to allow it to persist. 

If the applied field is increased from a value for which there are two 
solutions, then, the solutions approach each other until at some field 
value they coalesce. This point at which the solutions coalesce is a 
point of radial metastability since, at this point, the construction line is 
tangent to F. Thus, the solution is neither stable nor unstable to lowest 
order by the small variation argument of the preceding paragraphs. 
There are no solutions for applied fields greater than the field of radial 
metastability. Since in this case the inward forces dominate at all 
diameters, the domain collapses. (Since F everywhere lies below the 
asymptote d/h, the domains always collapse for applied fields greater 
than 4rM s .) In Fig. 2 where l/h = 0.3, the collapse point, Q, occurs 
at H/±ttM, = 0.283, d/h = 1.16. The sequence of the types of solu- 
tions which occur as H is varied depends only on the invariant signs 
of the slope and curvature of F. Therefore for a fixed value of l/h this 
sequence is independent of the value of l/h. From the slope and curva- 
ture properties of F the uniqueness of the collapse diameter and field 
for a given value of l/h may also be shown. However, the discussion of 
both this uniqueness and the detailed behavior of the domain diameter 
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as a function of applied field are conveniently postponed until general 
stability is discussed. 

The device implications of this subsection may be summarized by 
noting that cylindrical domains exist only in the presence of an applied 
field applied in a direction tending to collapse the domain and having 
a magnitude less than 4:irM e . 

2.2 General Stability 

The sign of the second variation of the total domain energy pro- 
duced by a weak variation in shape characterizes the stability of a 
cylindrical domain. Since only domains in equilibrium are of interest 
here, attention is restricted to variations of such domains from their 
equilibrium size and strictly circular shape. Since an arbitrary weak 
variation is describable by the expansion (2) with condition (3), 
the stability problem is reduced to the study of the coefficients of the 
terms in the energy expansion which are quadratic in the Ar n and A0 n . 
These coefficients are defined (up to some constant factor) as the 
stiffness matrix elements of the system, and these stiffness matrix ele- 
ments may clearly be classified as either radial, mixed, or angular stiff- 
ness matrix elements. 

In the equilibrium energy expansion (7) the only nonzero quadratic 
coefficients are the coefficients of the (Ar„) 2 , n ¥• 1. As required by the 
cylindrical symmetry of the system, the domain is completely metasta- 
ble with respect to angle. This is indicated formally in (7) by the 
absence of any nonzero terms in A0 n A0„, or Ar„A0 M1 . (The angular and 
mixed stiffness submatrices are identically zero.) Since no terms in 
Ar n Ar,„ for m ¥= n appear in equation (7), (the radial stiffness matrix 
is diagonal) the variation amplitudes are quasi-normal-modes of the 
system. Nonzero terms of order three and higher in the variation ampli- 
tudes and angles prevent the variation amplitudes being true normal 
modes. Finally, since the energy of a domain in an infinite plate is 
independent of the domain position, the translational stiffness of the 
domain is zero and therefore since the variation An corresponds to 
lowest order to a translation (see Ref. 1 Sec. 4.2.2), the coefficient of 
(An) 2 in equation (7) is zero. 

Since the second-order energy variation with respect to an arbitrary 
small amplitude weak variation is simply the sum of the energy 
variations from each normal component of the variation, the study of 
the stability of circular cylindrical domains reduces to the study of 
the stability of the domains with respect to size, Ar , and shape, Ar 2 
to Ar M . The sign of the corresponding stiffness matrix element deter- 
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mines the stability or instability of the domain with respect to the 
variation of a particular Ar„. From equation (7), the nonzero normal 
stiffness matrix elements are 



h 2 



4(2ttM:) W \ dr\ 
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drl ) 
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1 - <i 
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n ^ 2. (15b) 



The quantities in square brackets are termed stability coefficients. 
The domain is stable with respect to an arbitrary variation in shape 
when all of the stiffness matrix elements are positive. From equation 
(15), this occurs when the n = coefficient is negative and all the 
other stability coefficients are positive, 



and 



or 



[l/h - S (d/h)] < 



[l/h - S n (d/h)] > 0, 



n > 2. 



S (d/h) > l/h > S 2 (d/h). 
Since the stability functions have the property, 

S n+x (d/h) < S n (d/h), 



(16a) 



(16b) 



(17) 



(18) 



(at least up to n = 10 see Ref. 1) the condition for total stability 
reduces to 



S Q (d/h) > l/h > S 2 (d/h). 



(19) 



Since domain stability with respect to Ar and Ar 2 assures total stability, 
attention is largely restricted to the n = and n = 1 coefficients. The 
variations Ar and Ar 2 are termed radial and elliptical variations 
respectively. 

The numerical value of the stability coefficients and in particular 
the conditions under which the stability coefficients change sign, is 
determined by a graphical construction, an example of which is now 
given. The construction of the radial and elliptical stability coefficients 
for the case l/h = 0.3 is shown in Fig. 2. The elevation of the horizontal 
line W(R' represents the value of l/h. In order to determine the stability 
coefficients it is necessary to specify an operating diameter which, in 
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turn, is determined by the applied field. When the domain diameter is 
greater than the diameter at the intersection of the horizontal line 
and S 2 at (R' (applied fields less than that represented by the slope of 
the line W(R) the elliptical stability coefficient is negative and the do- 
main is unstable with respect to Ar 2 . The domain diameter at the 
intersection of the horizontal line and S 2 is called the diameter of 
elliptical metastability and the corresponding field (the field represented 
by the slope of the line W(R) is termed the field of elliptical metastability. 
When the diameter is decreased below the diameter of elliptical metasta- 
bility (by increasing the field above the field of elliptical metastability) 
the domain becomes stable with respect to elliptical deformation. When 
the domain diameter is decreased to d = 2h (the corresponding field 
being represented by the slope of the line W8) the magnitude of the 
radial stability coefficient (corresponding to the length of §>'§>'„) is approx- 
imately equal to the magnitude of the elliptical stability coefficient 
(corresponding to the length of §>'§>' 2 ). 

Increasing the field to the value corresponding to the slope of the 
line We, decreases the domain diameter further to the value at the 
intersection of the horizontal construction line with S n at C', so that 
the radial stability coefficient is zero and the force equation construc- 
tion line is tangent to the force curve. Increasing the field above the 
field value corresponding to the slope of We decreases the diameter 
even further so that the radial stability coefficient becomes positive, 
the radial stiffness matrix element becomes negative, there are no 
solutions to the force equation and the domain collapses. 

Since the construction line for the force equation is tangent to the 
force curve at the diameter of radial metastability, it is possible to 
construct the S curve from the F curve by plotting the loci of points 
whose ordinate is l/h and whose abscissa is the diameter at which 
the construction line is tangent to the F curve. Figure 3 illustrates 
four points of such a construction. Conversely, given the initial slope 
of F, the S curve may be used to construct F. 

In Appendix A of Ref. 1 (see also Fig. 2 here) it was shown that 
the stability functions are in general monotonic increasing functions of 
the normalized domain diameter, d/h, and monotonic decreasing 
functions of the radial periodicity, n (at least up to n = 10). From 
these properties of the S„ functions and the properties of the solution 
to the force equation discussed in Section 2.1, the following general 
stability properties may be deduced for any fixed value of l/h: When 
there is no applied field, the domain diameter is infinite, and the 
domain is unstable with respect to all variations for which n i? 2. 
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Fig. 3 — Construction of the S n (d/h) function as a sequence of collapse diameter 
solutions. 



When the applied field is small, the diameter is finite and the domain 
is stable with respect to all variation for which n is equal or greater 
than some n,„, while remaining unstable with respect to variations of 
r„ for which 2 £s n < n m . When the applied field is between the values 
of elliptical and radial metastability, the domain is stable with respect 
to all variations. Finally, for applied fields greater than the field of 
radial metastability, the domain collapses. (As noted previously the 
domain collapses for applied fields greater than the fields of radial 
metastability for any value of the domain diameter.) 

Note that except for Si each S„ forms the boundary between the 
regions of stability and instability with respect to the corresponding 
Ar n and that metastability with respect to each Ar„ (n ¥= 1) occurs only 
along the boundary between the stable and unstable regions. 

Even though the radial and elliptical stability functions bound the 
region of total domain stability, the threefold and fourfold stability 
functions, S 3 and S 4 , lie quite close to S 2 and it might be expected 
that when the bias field on a cylindrical domain is reduced in the 
presence of a small stabilizing coercivity the domain might run out into 
an initially three or fourfold figure. Such runouts have indeed been 
observed. 8 

When the bias field on an initially circular stable domain, repre- 
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sented for example by the line *WS in Fig. 2, is decreased, r b (d) will 
increase uniformly for all maintaining a circular shape until the point 
(R is reached. Increasing the domain radius beyond this value causes 
the domain to become unstable with respect to elliptical variations. 
Since for small variations, the variations in the expansion used here 
are normal modes, the breakaway to instability will not be coupled to 
any pure radial motion. Thus, when the breakaway occurs r b (d) will 
actually decrease along one half the length of the wall. Experiments in 
platelets which have just enough coercivity so that the initiation of 
the breakaway is observable confirm this somewhat surprising pre- 
diction. (See Ref. 4, pp. 1916-1917. The sequence of photographs on 
these pages were taken with an increasing bias field. The corresponding 
sequence for a decreasing bias field is similar.) 

2.3 Restrictions Placed on the Possible Region oj Device Operation by 
Stability Considerations 

Since the diameters of radial and elliptical metastability (and the 
corresponding applied fields) are the boundaries of the region of total 
stability, they are the boundaries of the region of possible device 
operation in the hard bubble mode. It is easily appreciated by inspec- 
tion of equation (15) that the domain will be metastable with respect 
to a particular &r n {n ¥= 1) if and only if the corresponding stability 
coefficient is zero, 

l/h - S n (d/h) = 0. (20) 

Since the S n are monotonic increasing functions of d/h, a diameter 
of metastability is uniquely defined for each value of n and charac- 
teristic length to thickness ratio, 

djh = S; 1 (l/h), ft »* 1. (21) 

The corresponding applied fields of metastability are then defined 
using the force equation (14) by 



5© 

4tM. d 



d n \_ F \h ) h] ' 



(22) 



While l/h is a normalized wall energy, it is also, obviously, the recipro- 
cal of the thickness measured in units of the characteristic length. 
Thus, the d n and H n are functions of the normalized thickness, h/l. 
In the remaining topics of this section it will be appropriate to think of 
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this parameter as the normalized thickness, since in a given material 
it is the thickness which is the accessible variable. 

In principle, the values of the d n and H n for each value of h/l could 
be obtained by graphical construction. Figure 2 may be taken as an 
example of such a construction for h/l = 3.33 and n = and 2. In 
this construction the intersection at G' represents <2 (0.3) = 1.16h, the 
intersection at (R' represents d 2 (Q.3) = 3.58/i, the slope of WC repre- 
sents H (0.3) = 0.283(4ttM.) and the slope of -W(R represents /f 2 (0.3) = 
0.196(47rM,). (The values quoted here were obtained numerically.) 

To better appreciate the way in which the requirements of static 
stability bound the possible region of device operation, several diam- 
eter and applied field functions of the normalized plate thickness 
have been plotted and their asymptotic forms in the limit of very 
thick or very thin plates have been computed. 

The values used in plotting the functions were obtained in the 
following way: The implicit equation for d and d 2 [equation (20)] 
was inverted numerically using the expressions for So(d/h) [equation 
(12)] and *S 2 (d/h) [equation (13)] to obtain d /h and d 2 /h as func- 
tions of h/l. (Parametric plotting may be used only for functions of a 
single d„.) The desired functions were then computed from d^/h, d 2 /h 
and the expressions for the H„ [equation (22)] and F(d/h) [equation 
(11)]. (See also Ref. 1, Sections A.l and A. 8.) The asymptotes of the 
various curves are obtained using the expansions for F(d/h), S (d/h) 
and S 2 {d/h) from Appendix B of Ref. 1, with the force equation 
where necessary (14) and the d„ defining equation (20). The small 
d/h expansions are carried to order (d/h) 2 and the large d/h expan- 
sions are carried only to the lowest order constant and log terms in 
order to facilitate the algebraic inversion of equation (20) . 

Since the validity of some of the assumptions of the theory increases 
with increasing d/h, each plot includes an arrow indicating the direc- 
tion of increasing dn/h along the plotted curve or curves. The marks 
at do/h = 1 and d->/h = 1 indicate the points at which the theory 
becomes definitely suspect. It has been found experimentally, however, 
that the theory does give reasonably consistent results for values of d/h 
which are somewhat less than one. 

Figure 4 is a plot of the diameters of radial and elliptical metasta- 
bility measured in units of the characteristic length, do/l = (h/l) (d /h) 
and d 2 /l = (h/l) (d 2 /h), as functions of the thickness measured in 
units of the characteristic length, h/l. These diameters have the 
asymptotic values 
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Fig. 4 — Diameters of radial and elliptical metastability, do and da , in units of 
the characteristic length, I, as a function of the thickness, h, in units of /. 



do/l re 0.412(fc/Q exp (rl/h), h/l « 1, 



and 



(23a) 
(23b) 
(24a) 



do/l wl.268(fc/0*i h/l»l, 
for the collapse diameter, and 

d 2 /l re 1.564(A/0 exp (irl/h), h/l « 1, 
and 

d 2 /l & 3.760(fc/0*i fc/J » 1, ( 24b ) 

for the elliptical runout diameter. The uppermost curve in Fig. 5 is a 
plot of the ratio of the diameter of elliptical metastability to the 
diameter of radial metastability, do/d = (c^//i)/(c?o//i) = 
(d 2 /l)/ (d /l), as a function of the h/l. This ratio has the aymptotes 



d 2 /d re 3.794, h/l « 1 



(25a) 
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and 



d 2 /d « 3.000, h/l » 1 (25b) 

which are easily obtainable from equations (23) and (24) . Figure 6 is 
a plot of the relative variation of the geometric mean of the diameters 
of radial and elliptical metastability with respect to a relative varia- 
tion in plate thickness dln| (c?„<£>)-|/dln|/i| as a function of h/l. 
This curve has the asymptotes 



dlnKdodtfl/dln 
and 



1.000 - S.U2(h/iy 



h/l»l, 



(26a) 



a In {(dod^l/d In |/i| « 0.500, h/l » 1. 



(26b) 



Inspection of these limiting diameter plots and asymptotic expres- 
sions yields the following: In a given material, the minimum stable 
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Fig. 6— Logarithmic derivative of the domain diameter with respect to plate 
thickness, d In d/d In \h\, for the bias condition d = (drfUi as a function of 
the plate thickness, h. measured in units of the characteristic length, I. 



diameter attainable for any plate thickness is ss3.9i at a thickness 
of ~3.3Z. The minimum diameter of elliptical instability is ^121 
attained at a thickness of ss4.2Z. For any given plate thickness, the 
range of diameters over which the domain is stable is small, the ratio 
of the diameter of elliptical metastability to the diameter of radial 
metastability being roughly equal to three for any plate thickness. 
The domain diameter is thus for all practical purposes determined 
once the plate thickness and characteristic length are known. The 
increase in domain diameter with increasing plate thickness in thick 
plates is quite mild being according to a square root law. The increase 
in domain diameter with decreasing plate thickness is, on the other 
hand, exponential for a thin plate. This variation is so rapid that the 
magnitude of the relative variation of the diameter of a centrally 
biased domain, d = (d d 2 ) i with respect to a relative variation in 
thickness increases according to the inverse of the thickness. 
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Figure 7 is a plot of the applied fields of radial and elliptical 
metastability, H /4:irM g = (h/d )[F(d /h) - l/h) and H 2 /^M 8 = 
(h/d 2 )[F(d 2 /h) - l/h], as functions of h/l. These curves have the 
asymptotes 



Ho/lirM. « 0.772 exp (-vl/h), h/l « 1, (27a) 



and 



# /4ttM. » 1.000 - 1.596(V0"*, h/l » 1, (27b) 

for the collapse field and 

H 2 /4ttM. PS 0.475 exp (—rl/h), h/l « 1, (28a) 

and 

H 2 /4tcM. & 1.000 - 2.660(fc/0~*, V* » !. ( 28b ) 

for the elliptical runout field. The two lowermost curves in Fig. 5 are 
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Fig. 7— Applied fields of radial and elliptical metastability, Ho and H 2 , in units 
of 4wM. as a function of thickness, h, measured in units of the characteristic 
length, I. 
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plots of the ratio of the collapse field to the elliptical runout field and 
the ratio of the difference of these fields to their average as functions 
of h/l. The asymptotes of these ratios are 

H /H 2 » 1.626, h/l « 1, (29a) 

and 

H /H 2 & 1.000 + 1.064(ft/0"*, V » li ( 29b ) 

for the field ratio and 

2(H - H 2 )/(H + H 2 ) » 0.477, fc/Z « 1, (30a) 

and 

2(# - ff a )/(#o + H 2 ) te 1.064(VO"*i M » 1, (30b) 
for the ratio of the difference of the fields to their average. Equations 
(29) and (30) are easily obtainable from equations (27) and (28) . 

Inspection of these limiting applied field plots and asymptotes 
yields the following: The field magnitudes decrease exponentially 
with decreasing plate thickness, rapidly becoming unmanageably 
small for very thin plates. The fields increase monotonically with 
increasing plate thickness toward the common value, 4wM 8 . The rela- 
tive variation of applied fields, allowed within the region of stable 
circular domains, is smaller than the allowed relative variation of 
diameters, approaching a constant for thin plates and zero for thick 
plates. 

Note that 4ttM 8 determinations are most accurately carried out in 
moderately thick plates (within the limits of the cylindrical wall 
approximation and the coercivity limits described in the next section) 
since here H /4ttM 8 (and H 2 /4:wM 8 ) is a weak function of I and is 
asymptotic to one. Measurements of the characteristic length, on the 
other hand, are best carried out in moderately thin plates (within 
coercivity limits) where the diameter is a strong function of I. The 
appreciation of the validity of these statements will be greatly en- 
hanced if the reader will try a few constructions on a plot of the 
F and <S„ functions. 

The implications of the foregoing for device design are summarized 
as follows: For a given value of h/l there is only an approximately 
three-to-one variation in the domain diameter permitted within the 
stable region so that given material and plate thickness, domain size 
is closely determined. The minimum domain diameter occurs for 
h/l ~ 4. For thinner and thinner plates, both the plate thickness 
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margins and the applied field magnitude decrease exponentially. While 
field margins are reasonable for h/l ~ 4 (H /H 2 = 1.41) [or 2(H — 
H 2 )/(H + H 2 ) ~ 0.34], they become vanishingly small in very thick 
plates. These considerations thus demonstrate that from stability con- 
siderations alone, for a given value of the characteristic length, there 
is both a preferred range of thicknesses and a preferred range of 
domain diameters (or applied fields) for device operation. The 
preferred range of each is centered about the preferred values 

K = 4Z (31) 

and 

d p = 81 (32) 

for which the bias field is 

H p = 0.279 (4ttM,). (33) 

These values will be shown to be preferred in another sense in the 
next section. 

2.4 Domain Diameter as a Function of Applied Field 

Since domain diameter is in practice only measurable for diameters 
between d Q and d 2 (applied fields between H Q and H 2 ), these points 
form natural endpoints for plotting the diameter as a function of 
applied field. Figure 8 is a plot of {d - d n )/(d 2 - d ), as a 
function of (H - H 2 )/(H () - H 2 ), for various values of the 
normal thickness in the infinite disk. The plots for finite nonzero 
values of the thickness were obtained numerically using the force 
equation (14), the defining relation for the d n [equation (21)] and 
the expressions for F [equation (11)], S [equation (12)] and S 2 
[equation (13) ]. In the limit of very thick plates, diameter and applied 
field are related by 

EaszH = ^A^doY h ^„ (34a) 

H — H 2 4 dd 

using the small d/h expansions of F, S and S 2 from Ref. 1, Appendix 
B, to order (d/h) 2 . In the limit of very thin plates, diameter and 
applied field are related by 



Ho - H 1 [ _ a\ | , 

H - H 2 ~ 1 - 1 exp (-!) L d V 



d 

do 



h->0 (34b) 
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Fig. 8— Relative diameter, (d — <£>)/((&. - do), as a function of relative applied 
field, {H — H 2 )/(H — H 2 ), in the infinite plate for various values of the plate 
thickness, h, measured in units of the characteristic length, I. 



using the constant and simple log terms of the expansions of F, S 
and S 2 from Appendix B of Ref . 1 . The curves for intermediate thick- 
nesses are seen to lie, in order, between the limiting curves and are all 
nearly linear except in the neighborhood of the collapse diameter. In 
the neighborhood of the collapse diameter, the slopes of all the curves 
are infinite. This is easily appreciated by expressing the derivative of 
the domain diameter with respect to applied field in terms of the 
radial stability coefficient. 

The derivative of the domain diameter with respect to the applied 
field is obtained by considering equation (14) to be continuously 
solved and then differentiating with respect to d to obtain 



4ttM 



1 (ddH . 1 „\ d „,,,,, . 



(35) 



Eliminating H with the force equation, using the equation for $ (12) 
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and rearranging yields 



4ttM. dd = d 1 



(36) 



So that indeed the derivative of diameter with respect to applied field 
approaches infinity as radial metastability is approached. 

2.5 The Effect of Wall Energy Anisotropy on the Size, Shape and Stability 
of Cylindrical Domains 

The domains observed in orthoferrites are never precisely circular 
but always have some degree of ellipticity. The present subsection 
relates this anisotropy in domain shape to more fundamental domain 
parameters. 

In the orthoferrites K u » 2irM] so that the magnetization lies rigidly 
along the plate normal and the wall width is narrow as compared to the 
domain diameter (1). The applied field energy, E H and the internal 
magnetostatic energy, E M , terms of the total energy expression (4) 
thus make no contribution to producing the domain anisotropy. There- 
fore the anisotropy results from an anisotropy in the wall energy 
density <r w . Considering again the wall energy density to be inde- 
pendent of wall curvature and the wall to be oriented with its normal 
perpendicular to the plate normal, i.e., a cylindrical wall, the wall 
energy density may be expanded as 

00 

<?w = &w + h 2 °"2n cos [2n(v — p 2n )] (37) 

n = l 

where v is the angle between the wall normal, N, and the x axis (see 
Fig. 9) and the expansion coefficients, a w , a 2n and v 2n are taken to be 
positive. The odd angular periodicity expansion coefficients are deleted 
from equation (37) since the energy of the system is invariant under 
time reversal while the direction of the wall normal (referred to the 
magnetization direction) changes sign. 

To describe the anisotropy observed in orthoferrites, it has proven 
sufficient to include only the average and two-fold terms in equation 
(37). In addition, if the plate is assumed oriented so that the wall 
energy is maximized when the wall normal lies along the x axis, the 
wall energy density is 

<T W = o w + \ A<T W cos 2v. (38) 

Although the method which now will be employed to calculate the 
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Fig. 9 — Coordinate system used in consideration of anisotropic wall energy. 

implications of equation (38) is clearly applicable to the more gen- 
eral energy expression (37), attention will be restricted to equa- 
tion (38). 

The total wall energy is 



E x 



-hf 



ds 



(39) 



where s is the arc length along the curve describing the domain shape 
in the plane. When the wall energy density (38) is substituted into 
equation (39) and the differential arc length, ds, and the wall normal 
orientation angle, v (see Fig. 9), are expressed in terms of 0, r h {6) and 
dr b (6)/d$, the total wall energy expression becomes 



+ v 



- 2^M^rcos20 -r 6 sin 20 j|?- 2 6 + \~) I > dO. 



dd \dd 



(40) 
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In the isotropic wall energy density case, <r w = a w , the first term in 
equation (40) is identical to the isotropic wall energy expression [equa- 
tion (8) of Ref. 1], the second term in equation (40) representing en- 
tirely the effect of the wall energy anisotropy. 

To evaluate the effect of the anisotropy term even to lowest order, 
it is necessary to obtain all the first and second derivatives with 
respect to the expansion coefficients, Ar„ and A0„ of the expansion of 
r b (6) (2) for the case of a strictly circular domain r b (6) = r . In the 
expressions for the derivatives, dAr„ and dA0„ are again abbreviated 
to dr n and dd„ and evaluation at r b (0) = r is denoted by a zero sub- 
script. The first derivative with respect to Ar n is 



BE W , f 2 \ 

— — = ha w I -] 
dr„ Jo 











•&)'. 




fin 




dr h 


r b 


dr n 


+ 


dd 


dr n 



[* + &)■] 



dd 



+ hh A<r„ 



f 



dn , dr b \de 

Tb dr n dd dr n 



[* + te)" 



-\ 



cos 20 



-2^(^cos20-r 6 sin20 
dr n \dd 



£ I I *■" O 



_ 2 dn d 



cos 20 -r 6 sin 20 ][ri + l-~ 



(am 



dd. (41) 



dd dr n _\dd 
Setting r b = r„ and dr b /d$ = and carrying out the integration yields 

fd^jA = 2vh _^ s to h A<rw 8n2 cos 2d 2 (42a) 

V dr„ U * 

where S mn is the Kronecker delta function. The remaining derivatives 
which are evaluated in an entirely similar manner are: 



as, 



de n / 



( d 2 E w \ 
\dr m drj 



dr m drj n r 



= 0, 



« 7 - 2 e 

= — ha w n 8 mn 



(42b) 



3tt 



+ — h Aa w mn{8 ni 6 ml cos 20, - S n . mi2 cos (n0 n — w0 m )}, 
4r 

(42c) 
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(J&rX - »' (42d) 



and 

(jf%3 = 37r/i A<r " *"' 5 " 2 Shl 29z ' (42e) 

The terms of equation (42) in d w are all identical to the corresponding 
terms in the isotropic case [equation (14) of Ref. 1] so that adding the 
terms of equation (42) in Aa w to equation (7) yields the total energy 
expansion for the anisotropic case. The domain is therefore in equi- 
librium when the force equation for the isotropic case (14) with <r w = a w 
is solved and the force tending to make the domain elliptical is some- 
how balanced. Inclusion of the effect of second order energy variations 
terms solves the latter part of the equilibrium problem to lowest order. 

Before proceeding with this solution, it is appropriate to comment 
upon the significance of the various energy terms. As required 
by translational symmetry [Ref. 1, Sec. 4.2.2] {dE 6iW /dr 2 ) = — 2?v 
(d 2 E AW /dr]) when X = 2 where E AW is the Aa w contribution to the 
energy. The (d 2 E w /dr 2 l ) term in the stiffness matrix thus has only 
kinematic significance. On the other hand, — d[Ar 2 (dE w /dr 2 ) ]/dd 2 = 
— Ar 2 (dE w /dd 2 dr 2 ) is a torque tending to turn an elliptical domain 
into the direction in which the force tending to make the domain 
elliptical is most positive. 

To solve the elliptical equilibrium problem, it is convenient to write 
the energy variation expression in the form 

AE T = -FX + |XSX + a (X) (43a) 
where 

X = [xi , x 2 , z 3 , Xi , x 5 , x a , x 7 , • • •] (43b) 

h [Ar , Ar, , A0, , Ar 2 , A0 2 , Ar 3 , Ad 3 , ■ • •] (43c) 

and where the elements of the force vector, F, and the stiffness matrix, 
S, are 



F 
and 



■ - -®. < 4M > 



- - \£tX ■ (43e) 
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Setting the gradient of the total energy in the space of the x, equal 
to zero and solving for the equilibrium yields displacements when third 
order energy terms are neglected 

X = S _, F (44) 

where S -1 is the inverse stiffness or compliance matrix. From the com- 
ments in the preceding two paragraphs and inspection of equations (43) 
and (44), it can be seen that stable equilibrium is obtained when: 
r is the stable solution to the isotropic force equation (14) with <r„ = a w , 
6 n = and A9 n = for all n (allowing now negative values of Ar„) and 
Ar n = for n = and n odd. The remaining even n Ar„ are determined 
from equation (44) written with respect to these variables only. Be- 
cause of the special form of the resulting F and S, inverting only a 
finite submatrix of S yields X to any finite order in Aa w . When this is 
done there results 

f - \ ft " *©]-(« + <T ) ■ 

*-ift-«ffln 



h S, \hIA ' / ' ' " u I. 



m+osm, ( 45b) 



f-<f)- ^^ (45c) 

where 

Al = AaJ^M 2 , . (46) 

Equation (45a) has also been obtained by E. Delia Torre and M. 
Dimyan. 9 

In the small Ar„ approximation of this calculation the wall anisotropy 
term has no effect on r or any of the odd n Ar n . The collapse diameter, 
do , and collapse field, H , are thus unchanged with respect to the iso- 
tropic case. At all diameters, however, the anisotropy pulls the initially 
circular domain out into an elliptical shape, the ellipticity being smallest 
at the collapse diameter and blowing up to infinity at the isotropic 
elliptical runout diameter (45a). At the isotropic runout diameter 
this calculation clearly cannot be applied to obtaining the reduction in 
the range of stability which wall anisotropy produces. It is also clear 
from the blow up of A>- 2 that any statement of such a reduction in 
the region of stability should be in terms of applied fields rather than 
domain diameters. In calculations of a w and Aa w from measurements 
of the major and minor axes of elliptical domains, it should be noted 
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that Ar 4 /r [equation (45b)] makes a small contribution to the length 
of each axis which is not removed by averaging. 

In platelets having the preferred thickness, h = 41, and biased to the 
collapse diameter, Ar 2 /r Q = 2M(Al/h) and Ar 4 /r = 4.65 (Al/h) 2 , while 
in the neighborhood of the isotropic elliptical runout diameter where 
Ar 2 -> oo, Ar 4 /r = 3.59(AZ//i) Ar 2 /r n . 

In rare earth orthoferrites, the anisotropy in wall energy is typically 
(within a factor of two) Acr w /a„ = 3% so that Al/h for h = 41 is only 
3/4%. The ellipticity at collapse is thus small. The contribution of 
Ar 4 is very small at all diameters (h = 41) and further, near collapse, 
the cylindrical wall and infinitesimal wall width approximations quite 
possibly produce larger errors. 

The implication of the foregoing for device applications is that wall 
energy anisotropy has no effect upon the collapse diameter and collapse 
field. As far as the elliptical runout conditions are concerned, careful 
measurements carried out on a low coercivity platelet of TmFe0 3 in 
which Aa w /a w = 3% show that the reduction in the ratio of the col- 
lapse to the elliptical runout field is only of the order of 1%. 10 

III. COERCIVITY AND MOBILITY 

The preceding sections treated only forces arising from reversible 
processes. The present section considers forces arising from the irre- 
versible processes describable by the wall motion coercivity and wall 
mobility. Several relations between wall parameters and domain 
parameters are obtained. In particular, this section computes the 
domain coercivity and mobility in terms of the wall coercivity and 
obtains the maximum allowable coercivity which permits device 
operation in the hard bubble mode. Inversely, it is shown how the 
wall mobility and coercivity may be obtained by domain measure- 
ments. 

3.1 Domain Dissipation 

The method used to take dissipative effects into account is to 
compute the power dissipation produced by a general variation in 
domain shape using the wall dissipation equation and then to set 
this equal to the power produced by the variation. By this procedure 
the equations for the various modes (translation, size change and 
deformation) are obtained. The dissipation equation approach is taken 
as a best first guess to the solution of coercivity problems. The reader 
is cautioned not to take any of the coercivity results too literally, since 
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coercivity is never uniform and in many cases depends on the direc- 
tion of wall motion. 11 

Wall motion in many materials is describable in terms of a wall 
motion coercivity and a wall motion mobility. 12-14 Such a description 
should be valid in any material for sufficiently low velocities and 
coercivities. Presently known uniaxial materials such as BaFe^Oio 
for which the velocity-drive field relation is not describable by 
coercivity and mobility at high drive fields 15 - 18 (they show roughly 
a limiting velocity) have low mobilities which reduce their usefulness 
in device applications. The present work therefore restricts attention 
to 180° domain walls in materials having the velocity-drive field 
relation 

_ f Ml0 (|tf L | - H c ), \H L \ > H c , (4?) 



0, \H L \ g H c 

where: v„ is the local wall velocity in a direction normal to the wall, 
H L is the total local field component parallel to the magnetization, and 
where equation (47) serves as the defining relation for the wall motion 
coercive field, H c , and the wall mobility fi w . The field H L includes 
the effect of all magnetic fields as well as any effective fields such as 
the "wall energy field," and may vary from point to point. The 
symbol, H r „ is used to distinguish it from H, used elsewhere to denote 
the spatially uniform z component of the applied field. For a planar 
wall in the absence of internal magnetostatic forces, H L is the applied 
field. Equation (47) is assumed to hold whether or not the wall is 
accelerating since for presently known materials wall inertia effects 
are negligible. 

Consider now a segment of 180° domain wall such as a portion of 
the domain wall shown in Fig. 1 and assume that H L is positive when 
it lies in the positive z direction. Under these conditions, the power 
input per unit area from the local field to an inward moving (decreas- 
ing r b ) segment of the domain wall is independent of the details of 
the magnetic configuration within the wall and is 

Pin = 2M, \H L v n \. (48) 

Since inertial effects have been assumed to be negligible, the power 
input must be equal to the power dissipated per unit wall area, p diB8 , 
so that eliminating H h from equations (47) and (48) results in 



p,ii 



= 2M,[h c \v n \ +— vly 



(49) 
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Integrating the dissipation density over the domain wall area yields 
the total dissipation 

P«„ = f 2M.\ H c \v n \ + — ojj Ida (50) 

•'wall L Mu> J 

where da is the differential wall area. For a circular domain, where 
i"b(0) — r o, this expression becomes 



= 2MMo f* [ff. Ivirl + — |v-i P | 2 ] 



dd (51) 



where i r is the unit vector in the radial direction. Now for a circular 
domain 

v .i = ^ 
V lr dt 



_ v l(^M\ drn , (dn\ dej 

& LA Sr n Jo dt "*" \ddj dt J 



= E cos [n(d - d n )] ft (52) 

where again dAr„ and dA6„ have been abbreviated to dr„ and dd n . 
Substituting equation (52) into equation (51) and carrying out the 
integration yields 

p -=^>-Hltl + fSltl + ™®] 

where NOi (dr„/dt) indicates the nonlinear coercivity coupling terms 
which appear even in lowest order. The nonlinear coupling terms 
tend in general to couple in additional modes even when only one 
drjdt is initially nonzero. An exception to this is the uniform radial 
mode of motion dr Q /dt which, because of its symmetry, may take place 
without coupling in the other modes. If the coercivity is negligible 
[H c <3C l/fLnidrn/dt)] then the Ar„ and A0„ remain uncoupled normal 
modes of the system. The consequences of the damping of the various 
individual modes will now be considered in the order n = 1, n = 0, 
n^ 2. 

3.2 Domain Coercivity and Mobility (n = 1) 

Consider an initially circular domain in which only d&r^dt is 
nonzero. In this case the component of the domain wall velocity 
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normal to the wall is 

v n = dr b (d)/dt = (d&rjdt) cos (0 - 0,). (54) 

Since the distribution of the wall velocity component normal to the 
domain wall for a circular domain propagating with velocity, v d , in 
the 0d direction is 

v n = i r v d = \v d \ cos (0 - d ), (55) 

a variation in Ar x may be identified with translation in the 8 d direction 
with 

|v„| = \dAr t /dt\. (56) 

Using equation (53), the power dissipated by a uniformly moving cir- 
cular domain is 

P di88 = *Mm\- H c \v d \+ - v°l (57) 

and, since P dl8fl = -v d -F d where F d is the drag force, the equivalent 
force for a domain moving with a nonzero velocity is 



= -TrMM \-H e + -\v d \\ 

Lit Hw J 



(58) 



where i„ = v d /|v d | is the unit vector in the direction of the motion. 
Notice that F d is an ordinary linear force which could be measured 
mechanically with the aid of a magnetic probe. 

Translational forces are most easily produced by gradients in the 
applied field. The power input to the domain from such a force is 
obtained by integrating the power input density, equation (48) , over 
the domain wall area. In the case of a uniform gradient the local field 
at the domain wall is 

B h = S L - i |AH| cos (0 - O ) (59) 

where AH is a vector orientated in the direction in which the bias field 
decreases most rapidly, 0„ , and has a magnitude equal to the max- 
imum difference in field across the domain (a gradient of magnitude 
\AH\/2r Q and direction 0„ + ir). In equation (59) H L includes: the 
bias field at the center of the domain, the "wall field" and the de- 
magnetizing field, all of which are independent of angle. All of the 
field components in equation (59), H L , H L and AH are to be understood 
as the z-averaged z components of the actual field since this is the 
quantity which interacts with a 180° cylindrical domain wall. With 
this understanding, integration over the wall area may be replaced by 
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multiplication by h and integration over 0. Assuming the domain to 
be exactly circular, the total power input to the domain is 

P ia = 2M.r h 

(•2' 

• / [H L - \ |AH| cos (0 - 6 B )]-[- \v d \ cos (0 - e d )] de 

Jo 

= Trr hM. AHv rf (60) 

which is just the expression for the power input to a dipole of strength 
(irr 2 h)(2M,) propagating in a gradient of magnitude |AH|/2?- . 

The domain will propagate in the direction in which the input power 
is greatest, (0„ = 9 d the direction in which the bias field decreases most 
rapidly), and the magnitude of the velocity will be such that the dissi- 
pated power (57) balances the input power (60) . Such a balance always 
occurs for v d = 0. When |AH| ^ (8/t)H c this is the only condition in 
which the balance is maintained. The domain velocity is therefore 



M = iU»l -\" 



I--HJ, \AK\>^H C , (61a) 

7T / T 



M = 0, |AH| ^ -H e . (61b) 



IT 



Comparison of equations (47) and (61) shows that it is possible to 
define a domain mobility and coercivity by taking AH as the driving 
field in terms of the wall mobility and coercivity as 

u. d = hn„ (62) 

and 

H cd = -H c (63) 

7T 

if AH is taken to be the drive field. 

Note that in the integral for the power input (60), the angle inde- 
pendent field term, ff L , does not contribute so that the effective power 
input density with respect to angle for a domain propagating down a 
bias field gradient is r hM. |AH| |v d | cos 2 (0 — 6 d ) which has the same 
angular factor as the mobility associated dissipation (50) and (55). 
The power input and power dissipated thus cancel (after z averaging) 
at each point on the perimeter of the domain for domains propagating 
under the influence of a uniform field gradient and viscous damping. 
In this case therefore there are no forces tending to distort the domain 
shape from circular. When coercivity is present (even perfectly uniform 



CYLINDRICAL MAGNETIC DOMAINS 



757 



coercivity), the power input and dissipation do not cancel locally and 
thus internal stresses tending to distort the domain from circular 
are present. 

When a bias field is applied to a circular domain which is not simply 
a uniform gradient, the affect of the additional nonuniformities may be 
taken into account by Fourier decomposition of the z-averaged z com- 
ponent of the applied field at the domain wall with respect to angle. 
When this is done the constant term determines domain size, the 
term translates the domain, and the nd terms, for n ^ 2, deform the 
domain. The procedure is, of course, only applicable if the domain 
shape remains near circular. 

Coercivity and mobility may be measured either by applying fields 
to walls 19 and using equation (47) or by applying field gradients to 
entire cylindrical domains 20 and using equation (61). In the platelets 
of material which have a coercivity which is so low as to make them 
useful for device applications it has been found convenient to use a 
second domain to provide the field gradient for the coercivity measure- 
ments. In this measurement two domains are brought together and 
then released, the coercivity being given by the formula 5 

tf c = (4ttM s ) |^ (64) 

where s is the center-to-center distance of the two circular domains. 
Equation (64) was obtained using equation (61) and the z component 
of the dipole field from the second domain at the center of the plate, 

2M.irr 2 h/s\ 

3.3 Domain Size in the Presence of Dissipative Processes (n = 0) 

The domain size variation mode Ar has two properties which are 
not common to the other modes: First, large variations from equi- 
librium may be considered using the force function, F(d/h), as well as 
small variations from equilibrium using the radial stability function, 
S (d/h). Second, arbitrary dissipation functions may be considered 
because the rotational symmetry of the motion removes any tendency 
for an initially circular domain to couple in other modes. 

3.3.1 Large Variations in Domain Size 

Setting all the clAr n /dt except dAr Q /dt equal to zero in equation (53) 
yields the dissipation produced by a domain size change 

P^-WU^+i^]^ (65) 
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where the upper sign is for an expanding domain. Setting the sum of 
the power dissipation and the rate of energy change from equation (7) 

f 2 = ft $ - «w**|J + t&- *©] d -f m 

equal to zero yields the differential equation for the domain diameter 
[dAr = dr n = id(d)] 

h d Ah) I , dH ± H c Jd\ , am 

= 1 + 1 „.„• -Ml • (67) 



2 Mu ,(4ttM'.) A d* fe ' h 4ttM. \ft 

The domain diameter thus relaxes toward a value which is a solution 
to the force equation (14) in which the bias field H has been replaced 
by a composite bias field, H ± H c . After this substitution is carried 
out, the equilibrium diameters are obtained as described in Section 2.1 
except that the sign of ±H C must now be determined in each case and 
now there is a small continuous range of stable solutions about both 
the stable and unstable zero coercivity solutions. There are again two 
solution diameters for each H ± H c when < H ± H c < H . In the 
present case, however, coercivity produces two stable ranges of solutions 
rather than one stable solution point and one unstable solution point. 
The large diameter solutions to the force equation for H + H e and 
H — H c bound the solution range which brackets the zero coercivity 
stable solution and similarly the small diameter solutions to the force 
equation for H -j- H c and H — H c bound the solution range which 
brackets the zero coercivity unstable solution. 

Figure 10 shows the graphical construction for the case l/h = 0.300, 
H/4-rrM, = 0.2544, H c /AttM, = 0.020. The zero coercivity stable and 
unstable solutions at d/h = 2.000 and d/h = 0.767 are marked S and "U 
respectively; the large diameter solutions for H + H c and H — H c at 
d/h = 1.527 and d/h = 2.468 are marked S + and S_ respectively and 
the small diameter solutions for H + H c and H — H c at d/h = 0.919 
and d/h = 0.685 are marked c \l + and 11. respectively. A domain having 
a diameter greater than that at S_ will relax toward the diameter at S_ 
as indicated by the arrow. Coercivity stabilizes domains having di- 
ameters between S+ and S_ . Domains having diameters between the 
diameter at S + and C U+ will relax toward S + as indicated by the arrow. 
Coercivity stabilizes domains having diameters between t U + and 11- . 
Small diameter coercivity stabilized solutions have been observed in 
the process of carrying out the mobility measurement described in the 
second paragraph which follows. 7 Domains having diameters smaller 
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Fig 10— Construction of solutions to the force equation in the presence of 
coercivity, l/h = 0.3, H/i*M. = 0.2544, H C /^M. = 0.020. 



than 11- will collapse as indicated by the arrow. The rate of change of 
domain diameter in any of the dynamic processes described is, of course, 
given by equation (67). 

It can be seen that the collapse diameter is independent of coercivity 
to the extent that the coercivity is uniform (this is only approximately 
true at best). A measurement of the collapse diameter, denoted by d m 
here, yields the characteristic length 



I = kSoidJh) 



(68) 



independent of the coercivity. If H c is then measured from the diameter 
hysteresis and this value is subtracted from the measured collapse 
field, the value of 4ttM. is obtained. 

Bobeck 17 ' 18 has developed a method for measuring mobilities which 
requires only the observation of static domain states. In this method the 
domain is biased to a stable diameter and then a short duration field 
pulse having a magnitude, H P , such that the total field amplitude 
exceeds the collapse field, H + H P > H , is applied. The pulse dura- 
tion is at first kept so short that at the end of the pulse, the domain 
has a diameter between the zero coercivity stable and unstable solu- 
tion so that it recovers its initial size. The pulse length is then gradually 
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increased until the diameter at the end of the pulse is less than the 
unstable solution diameter so that the domain collapses. 

Figure 11 is a construction illustrating the mobility method for 
l/h = 0.30, H/4:wM. = 0.254, H P /^irM, = 0.246 and zero coercivity. 
Under these conditions, (H + H p )/4tM, = 0.500 > H /4irM, = 0.283, 
and the solution diameters are d,/h = 2.00 > d /h = 1.16 > d u /h = 0.77 
where d, is the stable solution diameter and d u is the unstable solution 
diameter. In the figure the static collapse point is denoted by e. The 
domain is initially at the stable equilibrium point S. Application of a 
fast rise bias field pulse of amplitude H P takes the domain to S' where 
it collapses towards W. If the pulse is turned off at oY (above the un- 
stable solution point) ; the diameter increases until the point S is again 
reached. If, on the other hand, the pulse has sufficient duration to 
reach 2D' so that the domain state reaches 2), the domain collapses 
to W. In the insert, the field pulse shape used in making the measure- 
ment is shown as a function of time. 

When the domain velocity is proportional to the applied field and 




Fig. 11— Construction used in the pulsed bias field mobility method, l/h = 0.3, 
H/4rM, = 0.254, (H + H P )/4wM, = 0.500. 
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pulse rise time effects may be neglected, the mobility is 
h [ 

Mu> nrpA — T\/T J 



*d./h 



d A 



2T47tM, 



hu/h l_ . d H _ „(d 
h + h 4xM. \h 



(69) 



where T is the minimum pulse duration which results in the collapse 
of the domain. When coercivity must be taken into account the limits 
of the integral change as described previously. (At the beginning of the 
pulse the domain is at the point S + of Fig. 10.) If additionally nonlinear 
velocity-drive field relations and pulse rise times must be considered, 
it is not possible to obtain the mobility as a simple integral over the 
diameter and the integration is better carried out with respect to time. 

3.3.2 Small Variations in Domain Size from Equilibrium 
This section considers the effect of dissipative processes on domain 
size for the case where the energy variation expression (7) may be 
considered an expansion of the energy about the domain diameter 
which is a solution to the force equation (14). In this case the rate of 
change of energy with respect to time using (7) is 



dEj. 
dt 



_ dE T dAr 
dAr dt 

- -Me,** - d \j - 4i 



Ar 



dAr 
dt 



(70) 



Again equating the rate of energy decrease to the power dissipation 
(65) yields for the linearized differential equation for radial motion 



1 



dAr 



i)fi - *©] 



p w (4*M t ) dt 

The domain radius thus relaxes towards 



Ar 



T- 



JL 



AirM. 



Afo 



H, 



4wM, 



d/h 



l/h - S (d/h) 



(71) 



(72) 



if the domain is stable or, in the case of an unstable domain, coercivity 
stabilizes the domain for departures in radius from equilibrium up 
to this same value. In either case, stable or unstable, the relaxation 
time [defined by the time factor exp( — t/r)] is 



d/h 



M„4rM. l/h - S (d/h) 



(73) 
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The factor d/h[l/h - Soid/h)]' 1 in equations (71), (72) and (73) 
is a measure of the radial compliance of the domain relative to co- 
ercivity or mobility and will be called the radial relative compliance 
function. The value of this function may be obtained as the inverse 
slope of a line constructed on a plot of the force and stability functions. 
Such a construction is shown in Fig. 12 for l/h = 0.30, d/h = 2.00, 
H/4:irM, = 0.2544. The numerical slope of the line V?(R drawn from l/h 
on the vertical axis to the point on S at d/h = 2.00 is 0.0883 so that 
d/h[l/h — So]' 1 = 11.33. Thus, in this case, for a domain to have a 
diameter denned to within ten percent, the coercivity must be 
H c < 0.01 (4ttM.). 

At the collapse diameter, d , the relative radial compliance is, of 
course, infinite while at the other end of the range of stability (the 
elliptical runout diameter) where in the present example d 2 /h = 0.358, 
the relative radial compliance has the value 9.57. The minimum 
value of the relative radial compliance for l/h = 0.30 is 9.52 occurring 
at a diameter of d/h = 3.22. This behavior of the radial relative 
compliance at l/h = 0.30 is typical for thicknesses near the preferred 
value of h/l = 4 as can be seen from Fig. 12. In all such cases, the 
minimum compliance is achieved at some diameter less than d?, the 
compliance being nearly constant from d = d 2 down to (cWs) 1 [for 
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Fig. 12 — Construction of the radial and elliptical relative compliance functions 
for l/h = 0.3, d/h = 2. 
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l/h = 0.30, (dod 2 ) i /h = 2.04] then increasing rapidly to infinity at d^. 
The constancy of the radial compliance over a large part of the stable 
range is related to the linearity of the diameter-field curves away from 
the collapse diameter [see (36) and Fig. 8]. 

In Section 3.3.1, the diameters bounding the coercivity stabilized 
solutions to the equilibrium problem were computed for the case 
l/h = 0.30, H/4*rM 8 = 0.2544 and {d u /h = 0.767 d„/h = 2.00 and 
H c /4ttM 8 = 0.020 using the force equation in the presence of coercivity. 
As an illustrative consistency check, these bounding diameters will 
now be computed using the linearized equation (72). The diameters 
bounding the solution region bracketing the zero coercivity unstable 
solution obtained from this equation are 



!-*<■* 



fr *fetB -*&)]"'} - oj »* and °- 869 - 



and the diameters bounding the solution region bracketing the zero 
coercivity stable solution region are 



* fr * &t g - *®n -"•-"»■ 



j-™* 



The corresponding diameters obtained directly from the force equation 
were d/h = 0.685, 0.919, 1.527, and 2.468. The excellent agreement 
of the diameters bounding the zero coercivity stable solution com- 
puted by the two methods is again related to the linearity of the 
diameter-field curves. Note that the agreement of splittings of the 
diameters bounding the zero coercivity unstable solution computed the 
two ways is also quite good. 

3.4 Domain Shape in the Presence of Dissipative Processes (n ^ 2) 

The power dissipation produced by the motion of a single mode is 
from equation (53) 



P diB8 = 4*MM 



*''-+si^]^' »* 2 - (74) 

where the upper sign is for positive dAr„/dt. The rate of change of 
domain energy under these conditions is from (7) 

dE T = dE T dAr n 
dt ~ dAr n dt 
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Setting the sum of the rate of change of the energy and the power 
dissipation equal to zero yields the differential equation for the 
relaxation of the variation of a single Ar„, 

1 d£ . _ (n * _ Ml - slfj] & ± iJL. (76) 

m.4tM, dt K '\d/\_h WJ r ttItM, 

The domain shape variation amplitude thus relaxes towards 



|ArJ H e 



1 d/h 



4xM. ir (n 2 - 1) [l/h - S n (d/h)] I ' n ~ 2, (?7) 

if the variation is stable or, in the case of an unstable domain, 
coercivity stabilizes the domain for variations in amplitude up to this 
same value. In either case, stable or unstable, the relaxation time 
(defined by the time factor exp( — t/r)) is 

r 1 d/h 



' ' n t MM a (n 2 - 1) [l/h - S n (d/h)] ' U ~ 2 ' (78) 

Whenever the coercivity is effectively zero so that the mobility 
characterizes all dissipative processes, the normal modes of the domain 
remain decoupled (within the small amplitude approximation). Each 
of these modes relaxes according to equations (71) or (76). When 
coercivity must be considered, these equations become rather crude 
approximations because even perfectly uniform coercivity introduces 
nonlinear mode coupling. The nonlinear mode coupling is especially 
noticeable at the end of the relaxation of a single mode. The reason for 
phrasing the discussion of the effects of dissipative processes in terms 
of dissipation equations was to account correctly for the effects of 
coercivity to lowest order without being required to examine the 
coupling of the modes or the origins of coercivity. The results obtained 
do provide a general picture of the dependence of the effect of coerciv- 
ity on the various domain parameters and in particular they provide 
a measure of the dependence of the stiffness of the domains on these 
parameters. 

Equations (77) and (78) show that the entire n dependence of the 
residual distortion of a domain recovering from a fluctuation of a 
single mode and the relaxation time with which it recovers is contained 
in the stiffness factor (n 2 — l)[l/h — S n (d/h)]. The domain stability 
condition, S ^ l/h ^ S 2 , and the inequality (18) imply that [l/h — S n ] 
is a monotonic increasing function of n. Therefore the residual A?\, and 
relaxation time both decrease slightly faster than 1/n 2 (see Ref. 1). 
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Thus, to characterize the effect of coercivity in limiting the attainment 
of stable movable cylindrical domains it is necessary only to consider the 
elliptical shape variation mode (n = 2) in addition to the size (n = 0) 
and translation (n = 1) modes discussed previously. 

The factor (<L/3v)(d/k)[l/h - S 2 (d/h)]~ 1 in equations (76), and (77) 
and (78) for n = 2 is a measure of the elliptical compliance of the do- 
main relative to coercivity of mobility and will be called the elliptical 
relative compliance function. The value of this function is proportional 
to the inverse slope of a line constructed on a plot of the force and 
stability functions. Figure 12 shows such a construction for l/h = 0.30, 
d/h = 2.00, H/4ttM. = 0.2544. The numerical slope of the line *W6 
drawn from l/h on the vertical axis to the point on *S 2 at d/h = 2.00 
is -0.068 so that (4/Sir)(d/h)[l/h - S 2 ] -1 = 6.25. Thus a coercivity 
less than that value which defines the domain diameter to within ten 
percent (H e < 0.01(4irAf,)) defines the ratio of the difference in the 
ellipse semiaxes to their average (2Ar 2 /r ) to 12 percent. 

At the elliptical runout diameter, d 2 , the relative elliptical com- 
pliance is, of course, infinite while at the other end of the range of 
stability (the collapse diameter) where in the present case d n /h = 1.16 
the relative elliptical compliance has the value 2.22. The figure shows 
that this behavior is true for any value of l/h. In general the minimum 
value of the relative elliptical compliance occurs at the collapse 
diameter, the compliance increasing from the minimum in a regular 
fashion to infinity at the elliptical runout diameter. This regular 
behavior is in contrast to the behavior of the relative radial com- 
pliance. 

3.5 Restrictions Placed on the Region of Device Operation by Consideration 
of Dissipative Effects. 

The relative compliance functions appearing as factors in equations 
(72), (73), (77) and (78) contain the dependence of both the normal- 
ized residual distortion Ar„/r , and the normalized relaxation time 
r(/t»4irM a )/ro , upon n, I, h, and d, or H. Although the discussion which 
follows is phrased in terms of the residual Ar„/r , it should be kept in 
mind that the same remarks apply to the relaxation times scaled to 
ro/(/*w4irM a ), the time required to propagate the domain one radius 
when the maximum field difference across the domain, A//, is 8ttM„ 
[equation (61)]. 

In preceding subsections, the following has been demonstrated: 
Except in the immediate neighborhood of the collapse diameter, the 
domain diameter and bias field are, within the range of stability, 
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approximately linearly related (Fig. 8). In the neighborhood of the 
plate thickness previously termed preferred, h p = 41 [equation (31)], 
the radial compliance, considered as a function of the bias field, is 
roughly constant from d = (dad*)* to d = d 2 , and at d = (d d 2 ) i , 
the elliptical compliance is one half the radial compliance so that 
d = (d d 2 ) i represents a reasonable bias condition. Since for h/l ~ 4 
or in general for any plate thickness, the diameter and bias field 
ranges are relatively narrow (Fig. 5) the radial and elliptical compli- 
ance of domains suitable for device application may be completely 
characterized with respect to plate thickness by plotting the relative 
compliance functions with respect to h/l for the bias condition d = 
(dod 2 )i. 
The relative radial compliance, 

Aro/Vo d/h 



HJIttM, S (d/h) - l/h ' 






\,IVU,J 


» 1.08.WO*, 




h/l» 1, 


(79b) 


« 3.783 exp (irl/h) 


> 


h/l« 1, 


(79c) 


and the relative elliptical compliance, 








Ar 2 /r 4 . d/h 




! 


(80a) 


HJAttM, St S 2 (d/h) - 


l/h 



« 1.382(/i/J) i , h/l»l, (80b) 

« 1.606 exp (rl/h), h/l « 1, (80c) 

are plotted in Fig. 13 as functions of h/l for the bias condition d = 
(dody)*. The function values and asymptotic forms were obtained by 
methods which were used in Section 2.3 to obtain the diameter and 
field functions. The feature which distinguishes the relative compliance 
functions from the diameter and field functions is that the diameter 
and field functions bound the region of domain stability whereas the 
relative compliance functions provide a measure of the magnitude of 
the stability of the domains within the stable region. 

The constructions of the radial and elliptical relative compliance 
functions shown in Fig. 12 and described in Sections 3.3.2 and 3.4 
for h/l = 3.33, d/h = 2.00 may be taken as approximate construction 
for the values of these functions at h/l = 3.33 since (d d 2 ) i /h = 2.04. 

The minimum value of the radial compliance [d = (do^.)*] is ~ 7.9 
occurring at a thickness of h/l ~ 10.3 and the minimum value of the 
elliptical compliance [d = (dodo)*] is ~ 5.9 occurring at a thickness 
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Fig. 13 — Radial and elliptical relative compliance functions as a function of the 
thickness, h, measured in units of the characteristic length, I, for the bias condi- 
tion d = {d»dn)\. 

of h/l ~ 6.1. The compliance minima occur at somewhat greater plate 
thicknesses than the diameter minima. However, since the functions 
are quite flat bottomed, they increase only slightly in value from their 
minima to the thickness value previously termed preferred, h/l = 4.0 
[equation (31)]. {At h/l = 4.0, (d d. 2 )t/l ss 6.8. However, the pre- 
ferred values were taken to be h p = 41 [equation (31)], d p = SI [equa- 
tion (32) ] so that dp/hp = 2.0. The exact preferred values depend on 
the device structure in which the domain is located.} 

As in the example of Sections 3.3.2 and 3.4, the coercivity require- 
ment 

H c < 0.01 (4tM.) (81) 

will insure that a domain having the preferred thickness and diameter 
values will show radius function variations from the equilibrium 
radius of no more than ten percent. Just below h/l = 4.0, the increase 
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in the compliance functions is seen to become exponential so that the 
observation of stable moveable cylindrical domains in reasonably thin 
plates is experimentally unlikely. Since the relative compliance func- 
tions increase slowly with thickness for thick plates, cylindrical 
domains may be observed well into the region in which the assump- 
tion of cylindrical walls becomes dubious (see Ref. 1, Section 6.1). 

Although domains may be observed in thick plates, other considera- 
tions cause the realizable bit rate to decrease at least inversely with 
increasing thickness. These effects are as follows: The AH appearing 
in the domain velocity expression (61) is twice the n = 1 Fourier com- 
ponent of the 2-averaged z component of the applied field. Since the 
field gradient is applied from the surface of the plate and obeys La- 
place's equation, the variation in the applied field intensity decreases 
exponentially into the plate from its value at the surface. Thus for a 
given field variation intensity at the plate surface, the z-averaged z 
component of the intensity variation decreases inversely with increasing 
h/d for h/d > 1. Now for d = (d da)*, d/h = 1 occurs at h/l « 8 (see 
Figs. 2 or 13) and above this value d/h « 2.17(/i/0~ i [see Fig. 4 and 
equations (23b) and (24b)] so that for thicknesses greater than h/l « 8, 
the mobility with respect to the field gradient at the surface of the 
plate decreases according to (h/l)~*. Since in a given material the abso- 
lute domain diameter increases with increasing plate thickness in this 
region according to d/l « 2.17(fc/0*> the bit-rate decreases according 
to (h/l)" 1 . If additionally the maximum field difference at the surface 
of the plate is assumed to be some fraction of the difference of the 
collapse and runout fields, H — H 2 , then it is seen from equations 
(27b), (28b) and (30b) that the bit-rate decreases according to (h/l)'*. 

Thus in summary, consideration of the effects of dissipative 
processes even more strongly defines the neighborhood of h p /l = 4, 
dp/l = 8 as the preferred region than did considerations of stability 
only and additionally yields the requirement H c /4ttM 8 < 0.01 for the 
attainment of stable movable domains. 

IV. DETERMINATION OF MATERIAL PARAMETERS FROM PREFERRED DEVICE 
PARAMETERS* 

The preceding sections have provided preferred values of the plate 
thickness (31) and domain diameter (32) or bias field, and the least 



* Reference 21 includes part of the material of this section in a discussion of the 
relation of the M , and K» values of materials (available at that time) to the pre- 
ferred values of these parameters. 
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permissible value of the anisotropy constant (1). Additionally it was 
shown that the wall mobility acts to form the scale factor for time. 
For device construction it is desirable to specify the domain diameter, 
dp, from considerations of bit density and the resolution of mask- 
making and etching procedures while maximizing the bit rate. Once the 
domain diameter is specified, the desired characteristic material length 
is determined by equation (32) as I ~ d P /8, the thickness is deter- 
mined by equation (31) as h p ~ $d p , and the applied field is specified 
by equation (14) as H ~ 0.28(47rilf.,). It will now be shown that by 
adding assumptions about mobility and room temperature operation 
to conditions (1) and (31), it is possible to specify uniquely the three 
parameters A, K u and M s appearing in the energy density expression 
for the simplest uniaxial material 22 



Pe = A 



S+*"S 



+ tf u sin 2 d - H L M. (82) 



In equation (82), A is the isotropic exchange constant, is the polar 
angle (the angle between the magnetization and the z axis), <£ is the 
azimuthal angle, s is the distance through the wall, K u is the uniaxial 
anisotropy constant, M is the magnetization vector (|M| = M,) 
and H L is the sum of the applied and demagnetizing fields. Only Bloch 
walls in this simplest uniaxial material will be considered. Achievement 
of the coercivity condition (81) is a function of both intrinsic material 
properties and processing and will not be considered here. 

In Section VI of Ref. 1, it was shown that from several standpoints 

q m K u /2tM] , (83) 

the dimensionless ratio of the uniaxial anisotropy constant to the 
energy density of a volume containing a magnetic field of strength 
4nrM 8 must at least be greater than one, equation (1). From the 
definitions of q, I, and d p and the expressions for the Bloch wall width 
and energy in the simplest uniaxial material (82), l w = ir{A/K u )* and 
<t w = 4(AK„) § , 22 the ratio of preferred domain diameter to the wall 
width is 

f-^t. (84) 

If q were much less than one then the domain wall width would be 
larger than the domain diameter and clearly no domain of the type 
which has been considered here could exist. On the other hand, if q is 
very large the domain wall is very narrow with respect to the domain 
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diameter. For a given rate of flipping over of the spin systems com- 
prising the wall, the bit-rate will thus be inversely proportional to q. 
It may thus be expected that everything else being equal (it is not 
clear at the present time just exactly what is to be held constant) 
that materials with high q values will have low bit-rates in devices 
whose speed is drive field limited. A q value of approximately three 
may thus be termed preferred. 

The range of values of the exchange constant, A, occurring in ma- 
terials which may be considered for room temperature device applica- 
tions is quite limited as can be seen from the following argument: 
For a given structure and density of spin systems, the exchange 
constant is expected to be proportional to the Curie or Neel tempera- 
ture, T„. 23 For room temperature device applications, as a practical 
matter, T n must be above approximately 400°K. Since the highest ob- 
served T n are approximately 1000°K, the range of allowable T n and 
therefore the range of A in acceptable materials is nearly determined. 

F. B. Hagedorn, D. H. Smith, and F. C. Rossol have combined 
domain measurements of I with magnetometer measurements of K u 
and 4nM 8 to obtain the exchange constant in two materials. 24 They 
find for &m x Tb 1 - x Fe0 3 (x ss 0.55, T n = 661°K), A = 0.4 X 10" 6 
ergs per cm and for PbFe 12 - x Al x O w (x « 4.0, T n = 508°K) A = 
0.1 X 10 -6 ergs per cm, values which are apparently typical for high 
T n iron oxides. 

Since the exchange constant, A, is to be considered fixed, q has the 
preferred value three and d has the preferred value d p , it is appropriate 
to solve for the magnetization and anisotropy constant in terms of 
these quantities, 

4ttM. = S2(2irqA) i /d v (85) 

and 

K u = 256 Aq 2 /d 2 p . (86) 

If A = 4.0 X 10~ 7 ergs per cm, q = 3 and d p = 10~ 3 cm (approxi- 
mately one mil bit spacing) then 4tM. tt 80 Gauss and K u « 900 ergs 
per cc which are both numerically small. 

Maintaining the values of A and d p but considering q as variable (85) 
becomes 4:irM , « 50 \/~q. Thus q any within two orders of magnitude 
of the preferred value produces a value of 4:irM, of one kilogauss or 
smaller. Since the magnetic moment per spin system and the volume 
of the individual spin systems are approximately constants and since 
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the exchange interaction is of short range, low saturation magnetiza- 
tions cannot be achieved by dilution of the spin systems with magneti- 
cally unordered materials. Spin system density thus must be high 
while the net magnetic moment per unit volume remains low. Some 
sort of antiferromagnetic character is thus required in magnetic ma- 
terials which are candidates for use in cylindrical domain devices. 
Typical systems are ferrimagnets, and canted antiferromagnets. 

The K u value of 900 ergs per cc is quite low for a symmetry allowed 
intrinsic uniaxial anisotropy constant. 21 When low anisotropics are 
obtained by operating near the reorientation temperature in canted 
systems 25 or near the Neel temperature, the material parameters tend 
to be undesirably temperature dependent. Since uniaxial anisotropy 
energy densities of the required value may be induced, it appears 
that the use of materials having induced anisotropy, such as the 
recently announced garnets, 17,18,26 appears quite promising. 

Having determined that A is to be considered fixed and that K u 
and M a have preferred values, this section concludes by showing the 
dependence of several overall device parameters on these parameters 
and the mobility. The device parameters considered are the bit 
density, bit rate, the domain flux which is important in Hall effect 
detectors 27 and the induced voltage which is important in wire pickup 
loop detectors. 

In typical domain devices, the bit positions form a square array 
with the bit spacing being three to four domain diameters (see figures 
of References 4 and 5). If the bit spacing is assumed to be Zd p , then 
the number of bit locations per square centimeter is 

Pb = (3d p )" 2 = 6.9 X lO-'ilTrM.y/AK,. . (87) 

Since the bit density is such a strong function of 4ttM 8 the magnetiza- 
tion is nearly determined once a bit density is specified. 

The difference of the collapse and runout bias field for a plate of 
the preferred thickness h = 4J is H — H 2 ~ 0.1 (4irAf 8 ). Retaining 
the assumption of a bit spacing of Sd p and assuming that the device 
structure continuously maintains a field difference across the domain 
of 0.1(47rM g ) the bit rate in bits per second is 

j b = v b /M v = M .47rM./60d p = 4.1 X lO^n^irMyiAKX 11 , (88) 

linear in the mobility and again dominated by the 4rrAf g dependence 
of dp. If it is assumed that a pickup loop intercepts one half of the 
flux emerging from the magnetic charges forming the upper surface 
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of the cylindrical domain, then the flux change produced by moving 
a cylindrical domain under the loop in gauss square centimeters is 

* - lfi(4rM.) = 6.4 X 10 4 AK^M.y 3 . (89) 

o 

If it is assumed that this flux change takes place in the time required 
for a domain driven by a field difference of 0.1 (4wM a ) to propagate 
a distance d p , the induced (MKS practical) voltage is 

V = 7^ m.(4tM.)" d v X lO" 8 = 7.9 X \QT\MK$ . (90) 
lbU 
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